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1. Pointing Accuracy of the Low-Cost Large Space Telescope Due to 
Noise and Quantization 

1-1 . Introduction 

The objective of this study is to conduct an investigation on the 
pointing stability of the low-cost Large Space Telescope (LST) system. 

The low-cost LST is characterized by the use of reaction wheels for 
the generation of control torgues. Because of the critical reguirement 
on the pointing accuracy of the LST5 the nonlinear frictional characteristics 
of tte bearings of the reaction wheels cannot be neglected. It is well 
known that the nonlinear friction can cause limit cycles in a closed-loop 
system. 

Another possible source of pointing error in the LST is due to the 
effect of quantization and sensor noise. Since the LST is a digital 
system, D-A and A-D converters, and sensors for positional and rate 
feedbacks are used. Sensor noise and amplitude quantization will also 
cause pointing error in the LST. In addition, quantization is a 
nonlinear phenomenon so that it may also cause self-sustained oscillations 
in the closed-loop system. 

The dynamic modeling of the single-axis LST is described in this 
chapter. Several methods of evaluating the attitude error of the digital 
LST due to quantization and noise inputs are given. 
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1 - 2 . Dynamic Models of the Low-Co_st LST System 

The dynamic model of the single- axis low-cost LST system with 
sampled data is shown in the block diagram of Fig. 1-1. The rigid 
body is represented by the double-integrator transfer function. The 
controller is formed by proportional, rate, and integral feedbacks of 
the vehicle attitude. The nonlinear element N in the reaction wheel 
dynamics represents the rolling friction, and its functional description 

is given by the well-known Dahl model. 

The definitions of variables and the values of parameters and 

constants are tabulated In Table l-l [11 • 
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Table 1— !• 

Reference input command 
Body attitude of LST 
Body rate of LST 
Torque command of reaction wheel 
Torque output of reaction wheel 
Angular velocity of reaction wheel 
Angular displacement of reaction wheel 
Frictional torque of reaction wheel 
Attitude error 

Proportional gain of controller 
Rate gain of controller 
Integral gain of controller 
Feedback resistance 
Voltage amplifier gain 


1.65 X 10^ N-m/rad 
3.71 X 10^ N-m/rad/ sec 
7.33 X 10^ N-m/ rad/sec 
0.484 ohms 


10000 





Figure 1-1. Digital 1ow«cost 1ST lystM (slnglf aitft). 
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m 


Motor torque constant 
Back emf constant 

Moment of inertia of reaction wheel 

Moment of inertia of vehicle about 
pitch axis 

Motor time constant 

Motor resistance 


0.484 N-m/amp 
0.484 volt/rad/sec 
0.2 Kg-m^ 

41822 Kg-m^ 

0.002 sec 
10 ohms 


Although the low-cost LSI system has digital control, it is 
informative to analyze the system of Fig. 1-1 first without the sample- 
and-hald device. Figure 1-2 shows the signal flow graph of the continootis- 
data LST s^stea. The characteristic e<^ti<m of the systen is detemined 
from Fig. 1-2, 


4 - 1 ■» RpGjG} ♦ KgK^jSj 




^ 8r-6i6f,G.iG*G| 


'2”3"4 


F«1»2«3«4«5 


( 1 - 1 ) 


where 


G, = Kp t K,/s 

«2 = V* 

r = 1 

3 R + (1 + TjjjS) 



G 
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6 
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1 
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and N denotes the analog describing function of the reaction wheel nonlinearity. 



Figura 1-2. Stgi 
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For the linear model of the 1ST system, we set N = 0, and the 
characteristic equation becomes 

X J R-pS^ + j RtS^ + )S^ + RpKpK^s^ 

m V T vT 'vFA j-ka 

. K^RpKpS K^KjRp = 0 


(1-2) 


It is of interest to investigate the dynamics of the reaction wheel 
and the vehicle. The open-loop transfer function between and is 


<l> 


RcI^a 

r'-l T-^ ^ 

^ V ^ V * 


(1-3) 


Substituting the system paraiiwters into Eq. (1-3) gives 

h ^ 5.7865 

s^(s^ + 500s + 242058.5) 

5.7865 (1-4) 

s^(s + 250 + j423.74)(s + 250 + J423.74) 

Thus, the reaction wheel is shown to have relatively fast dynamics. 

Substitution of the values of the system parameters into tq. (1-2), 
and simplifying, the characteristic equation of the linear 1ST system is 
written 

s^ + 500s^ + 242058. 5s^ 2.1468 x lO^s^ + 9.5476 x lO^s 

+ 4.24145 X 10® = 0 (1-5) 


The roots of the characteristic equation are: 
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s = -0.49659 

s = -4.22743 + j4. 25123 

s = -4.22743 - j4. 25123 

s = -245.524 + j421 .118 

s = -245.524 - J421.118 

Note that the damping ratio of the dominant complex roots is 
0.705, and the natural undamped frequency is 6 rad/sec or 0.954 Hz. 

These parameters are achieved by selecting the controller constants K^, 

Kp and Kj at the indicated values. However, the poles of the reaction 
wheel dynamics at s = -250 + j423.74 and s = -250 - J423.74 are only 
slightly affected by the body controller and they account for the charac- 
teristic roots at s * -245.524 + j421.118 and s * -245.524 - j421.118 of 
the overall closed-loop system. Since these fast roots are very far away 
from the dominant ones, this means that for all practical purposes the 
dynamics of the reaction wheel can be neglected as far as the linear 
system is concerned. Figure 1-3 shows the block diagram of the simplified 
continuous-data low-cost LST system, and the digital system is shown in 
Fig. 1-4. 

The closed-loop transfer function of the continuous-data system of 
Fig. 1-3 is 








(1-6) 


For the digital system of Fig. 1-4, the closed-loop transfer function 


is written 
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$b(z) 

^^(z) 1 + 6^,Iz) + G^(z)Gg(z) 


where 



(s) 



0-7) 


-1 




KjS + Kj 


K z + K J - K 

= _£ I E. 

z - 1 


0 - 8 ) 


-1 


= 0-2 )| 


r 1 1 


V' 


- rUvO 

2J^(z - \) 


7 


0-9) 


-1 




f 1 1 


\ ^ J 


h 


Jy Z - 1 


0 - 10 ) 


Equation (1-7) is simplified to 


T^(KpZ^ ■»• K^Tz •»• K^T - Kp) 




2J^z^ + (T^Kp + 2K^T - 6J^)z^ + (6J^ - 4KgT + T^KjT)z + (2K^T + KjT^ - 2J^ - KpT^) 


( 1 - 11 ) 


The characteristic equation of the system is 


2J^z^ + (T^K + 2KgT - 6J^)z^ + (6J^ - 4KgT + t\j)z 


+ 2KgT + KjT^ - 2Jy - KpT^ = 0 (1-12) 


Or, 
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_ CO 

83644z^ + (1.65 X loV + 7.42 < lO^T - 2.50932 x 10 )z 

+ (2.50932 X 10^ - 14.84 x lO^T + 7.33 x 10^T^)z 

+ (-8.3644 X 10^ + 7.42 x lO^T + 7.33 x loV - 1.65 x loV) = 0 

(1-13) 

The characteristic equation roots are tabulated below as functions of 
the sampling period T: 


T (msec) 
0.1 
0.5 
1 
5 

10 

25 

50 

100 

120 

150 

170 

180 

200 

220 

230 


Roots 

1, 0.998 ± j0.0025 
1, 0,998 ± j0.0022 
U 0.995 t #.0043 
0.998, 0.979 ± j0.0208 
0.995, 0.957 ± jO.0411 
0.987 , 0.^ t #.0988 
0.975, 0.766 ± jO.1825 
0.950, 0.483 ± jO.2901 
0.940, 0.356 jO.300 

0.925, 0.150 ± j0.260 
0.915, 0.0034 ± j0.144 
0.910, 0.049, -0.195 
0.900, 0.107, -0.571 
0.890, 0.1116, -0.9078 
0.885, 0.1085, -1.077 (unstable) 


T. 


The root loci of Eq. (1-13) are sketched in Fig. 1-5 as a function of 
For small sampling periods, the characteristic roots are all located near 
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the z = I point. The linear digital 1ST system becomes unstable when 
T exceeds approximately 225 msec. 



1-3. Effects of Quantization on Pointing Stability of the Low-Cost 1ST - 
Limit Cycle Conditions 
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Quantization occurs at at least three places in the LSI system. Two 
are at the displacement sensor and the rate sensor where A-to-D converters 
are used. A quantizer is also needed at the control torque input to the 
reaction wheel since D-to-A conversion is effected there. In addition, if 
the integral control Kj/s is implemented digitally, quantization should be 
considered in the digital controller as well. Figure 1-6 shows the digital 
LSI system with quantizers. The z-transfer functions are defined in 
Eqs. (1-8), (1-9), and (1-10). The quantizers at the displacement 
sensor, the rate sensor, and the reaction wheel control torque are denoted 
hy Q , Qd» Qt» respectively. The input-output relation of a quantizer 
is shown in Fig. 1-7. The quantization level is represented by h. 

Ue shall analyze the effects of quantization on pointing stability 
or accuracy of the LST by means of three different methods. The first 
method utilizes the deterministic approach and establishes a lease upper 
bound on the pointing error due to quantization. The second method relies 
on treating the quantizer as a noise source, and statistical analysis is 
applied. The third method is also a statistical approach which represents 
the quantizer by a linearized gain Kgq(z)- 

It should be pointed out that a system with quantization is a nonlinear 
system and its behavior cannot be predicted by linear theory. One of the 
well-known phenomena of a nonlinear system is that sustained oscillations 
may occur. When a digital system has several quantizers, it is extremely 
difficult to predict the condition of self-sustained oscillations. 


jf-1 !?r*5!rv«iff?s» 
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To illustrate the effects of quantization, and how quantization can 
cause sustained oscillations in an otherwise stable linear system, let 
us refer to the digital systems shown in Figs. i-8 and 1-9. The difference 
between the two systems in Figs. l-8a and l-9a is that the former has 
negative feedback and the latter has positive feedback; but both systems 

are stable. 

For r(k) = 0, both linear systems have zero steady-state values for 
c(k); that is, c(k) = 0 for k-x», for arbitrary initial state c(0). We 
shall show that when quantization is considered, the system in Fig. l-9b 
has a steady-state error, whereas the syst«i in Fig. l-8b exhibits a 
sustained oscillation. 

Let the quantization level h be 2, (Fig. 1-7). The state equation 
of the system in Fig. l-9b is 

c(k+l) = Q[0.9c(k)] 

For c(0) = 10, it can be easily shown that c(k) = 2 for k > 8. 

The state equation of the system in Fig. l-8b is 

c(k+l) = Q[-0.9c{k)] 

For c(0) = 10, c(k) = 2 for k > 8 even, and c(k) = -2 for k > 8 
odd. Thus, the state of the system oscillates between -2 and +2. 
























1 
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1-4. Effects of Quantization on Poi n ting Stability of the Digital Low^ 

Cost 1ST - Least Upper Bound Q u antization Error 

in this section the effects of quontizution on the low-cost LST 
system are investigated using a deterministic approach. The method 
of analysis is based on the "worst" error condition due to quantization. 

In general, the analysis gives a conservative estimation of the quantization 


error. 

Since the quantization error has a maximum bound of ±h/2, the "worst" 
error due to quantization in a digital system can be studied by replacing 
the quantizer in the state dij^ram by a branch with unity gain and an 
extenml iwise source wftii a sl^wil saignitude of ±h/2. The block dia^a* 
of the digital LST system with quantizers shown in Fig. 1-6 is redr»m 
in Fig. 1-10 with the noise sources. The transfer functions %Cz). 

and G (z) are defined Tn Eqs. (1-8). (T-9), and (1-10), respectively. 

c 

The z- transform of the body attitude of the LST due to the three 
quantizers when = 0 is 




hpi h^ hjv 

± -^G.(z) ± -y ± 

1 + G^(z)6g(z) + gJzT 


0 - 16 ) 


Substitution of Eqs. (1-8) through (1-10) into Eq. (1-16), and 
simplifying, we have 


z * (^->-X^(z-l)]T^(z-H) 

V^> ° (z.,)3 + t2(z. 1)(K z.K,T-K ) + 2TK„(z-l)‘ 


0-17) 


i 
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The steady-state error of the body attitude due to the quantization effects 
is obtained by applying the final-value theorem to Eq. (1-17), (if the 
system is stable). 

Thus, 

tiiT.4>a(kT) = ftm(l - z'^>fg(z) = 
k-x» 2^1 

It is interesting to explore the significance of this result on 
the error due to quantization. Firstly, the quantization error at the 
displacement sensor is propagated through the system without change in 
anpli tilde. Secondly, the errors due to the torcpie and rate sensor 
quantizers are completely eliminated at the output position. This is 
attributed to the integral control Kj/s in the forward path. 

Vigitat JmpZejne.nt£it^n fohvxvtd ContfUftttA 

If the proportional-plus-integral controller is implemented 
digitally, the transfer function G^(z) becomes 


Substituting G^(z) from Eq. (1-19) into Eq. (1-16), we can again show 
that the steady-state error in <I>g(kT) due to the three quantizers is 
±hg/2, and the errors due to and are completely eliminated. 

In reality, the digital implementation of the controller should 
also include quantization in the digital process. Figure 1-11 shows 
the block diagram of the LST system with quantizations also considered 

The z- transform of the body attitude of the 


in the digital controller. 
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LSI due to the five quantizers when = 0 is 




fif, h, 

CH C 


h 




I) 


z - 1 


1 + G^TYTGgTzt^G^ 


o 


Gb(z) 


(1-20) 


where G^(z) is given by Eq. (1-19) and Gg(z) and G(.(z) are given in 
Eqs. (1-9) and (1-10), respectively. 

Applying the final -value theorem to Eq. (1-20), we have 


^^4>g(kT) = 

k-+® 


+ h + 
■ ”2 " 


hgKiT _ 


2KjT 



0 - 21 ) 


It is interesting to note that the introduction of the integral 
control eliminates the noise signals that enter at all points after the 
integral control in the control loop; however, the digital implementation 
of this control in turn produces a quantization error which is bounded 
by ±h2/2KjT, where h2 is the quantization level. 




1 - 5 . Effects of Quantization on Poi n tinq Stability of the Digital, Lowj! 
Cost 1ST - Equivalent Noise So u rce. Statistical Analysis 
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In this section the pointing errors uf the 


4 4 ^ r i. - 1 

u 1 I ca I 


I CT 
L- J » 
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to quanti Zdti ons in the displacement, rate, and torque channels are 
investigated by statistical means. The rms (root-mean-square) error 
in due to quantizations that are represented by equivalent Gaussian 
noise sources with zero mean, or white noise sources, is determined by 
setting = 0. The results are then compared with those of the continuous- 
data 1ST obtained in reference [1]. 

The block diagram model of the continuous-data LST with the <^antizers 
replaced by equivalent noise sources is shown in Fig. 1-12. The equivalent 
digital syst©n is shown in Fig. 1-13. It is assuned that the equivalent 
noise sources that represent the quantization operations are white, so 
that their power spectral density functions in both the s and the z 
domains are constants. Therefore, 


$^(z) = <J>p = power spectral density of displacement quantizer 

" 'I’d ” power Spectral density of rate quantizer 

$^(s), $.^(z) = = power spectral density of rate quantizer 

Let the rms attitude error of the continuous-data LST due to <I>p 
be represented by Ogp. The rms errors are given by the following 
relations: 

Displacement : 


°BD 


f 

1 

2irj 


• -j-» 


Mjj(s)Mjj(-s)ds • 


m 


1/2 


( 1 - 22 ) 






Figure 1-12. Continuous-data 1ST with guanttzers repres«ited hy 
equivalent white noise sources. 



Figure 1-13. Digital LST with quantizers represented by equivalent 
white noise sources. 
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where Mp(s) denotes the transfer function between and 'I>g{s), 


4>n(s) 

" $^TsT 


K^s . Kj 




(1-23) 


Similarly, 

Rate : 

^BR ^ 


Wf" 

J-Joo j 


1/2 


(1-24) 


Torque : 


'BT 


j®® 

2^ I • *T 

J -Joo 

L / 


1/2 


0-25) 


where 


-Mg(s) 


4>r(s) 


Jvs’ * * V* * '‘l 


(1-26) 


This result is obtained with the assumption that the quantizer noise 
is not multiplied by the sensor gain Kj^, as in reference [1]; otherwise, 
the right-hand side of Eq. (1-26) should be multiplied by Kj^. 

The total rms attitude error due to all three quantizers is simply 
the sum of the errors due to each noise source acting alone; that is. 


'^B ^ °B0 * *^BR * °BT 


(1-27) 


From Eq. (1-26) it is easy to see that Og,^ and Ogy have the same magnitude. 
The line integral of the form of Eq. (1-22) can be evaluated by 


I 


2A 


contour integration, and then the residue theorem. An integral table is 
also available for the evaluation of the contour integral. Using Eqs. 
(1-23) and (1-26), the ms attitude errors of the continuous-data system 
are obtained as 


^BD 


r + K,K„ 


'BT 


2(K„Kp - K,0^) % 


1 


1/2 


1/2 


'BR 




= 0.927 X 10'®*/5T 


(2.5748cf>p)^'^^ = 1.605/<J>p 


-12. ^l/2 


= (0.8598 X 10"'^$.^) 

= 0.927 X 10"®y^ 


For the digital LST, the rms attitude errors are; 
Displacement : 

iV2 


7* = 

^BD 


gU|Ho(2)MD(z-')2‘''dz • 


Rate: 


^* = 
^BR 




1/2 


Torque : 


s 

’bt 




1/2 


(1-28) 


0 - 29 ) 

(1-30) 


(1-31) 


(1-32) 


(1-33) 


The transfer functions, Mq(z), M^(z), and M.^(z) are determined from 


Fig. 1-13. 
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MqCz) 


Ga(z)Gb(z: 


1 + G^(z)6g(z) + G^(z) 


(1-34) 


Gg(z) 

-M^(z) = V\j(z) = 1 + G^(z)Gg(z) + G^{z) 


(1-35) 


where 


60 ( 2 ) 


= r(z + 1) 

2J,(Z - 1)‘ 


* j,(z - t) 


(1-36) 


0 - 37 ) 


and 


K Z + K,T - K 

G (zl = -B I B. 

z - 1 


(sample-and-hold and 
analog controller) 


0 - 38 ) 


G,(z) 


(2Kp + TK^)z TKj - 2Kp 


2(z - ly 


(digital implementation (1-39) 
of controller) 


Substituting Eqs. (1-36) through (1-39) into Eqs. (1-34) and (1-35), 
we get the following transfer functions which are used for the computation 
of the rms attitute errors in Eqs. (1-31) through (1-33). 

With zero-order hold and analog controller, 


T^CKpZ^ + KjZ + (KjT - Kp)] 
^1 


Md(2) = 


(1-40) 


26 


-^^(z) = M^(z) 




(1-41) 


where 

(-6Jv 2KpT + t \) z ^ + (6J^ - 4Kj^T + KjT^)z 
+ (-2J^ + 2K^T + KjT^ - KpT^) (1-42) 

With digitally implemented controller. 


T^[(2K + TK,)z^ + 2TK.Z + TK, - 2K 
M,(z) » E 1 I 1 E 


( 1 - 43 ) 


= H,(z) = ’> 


( 1 - 44 ) 


where 

A 2 = 4J^z^ + (-12J^ + 2T^Kp + t\j + 4K^T)z^ 

+ (12Jy + 2T^Kj - 8K|^T)z + (4Kj^T + T^Kj - 2T^Kp - 4J^) 

(1-45) 

A tabulation method or a numerical method [3] can be used to 
evaluate the contour integrals of Eqs. (1-31), (1-32), and (1-33), once 
the values of the parameters of the transfer functions are known. Table 

1-2 gives the values of (<J>Bnns^D‘ ^'^BnnsH ^‘*’Brms^R various 
values of the sampling period T. 


1 

I 
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Table 1-2. 

T Sample-and-Hold with Digitally Implemented 

(msec) Analog Controller Controller 



^BD 

^BT 

°BD 

^BT 

0.01 

0.0052 

2.932 X 10"^ 

0.0052 

2.932 X 10“^ 

0.1 

0.0161 

9.274 X 10'^ 

0.0161 

9.274 X 10'^ 

1.0 

0.0508 

2.936 X 10"^ 

0.0508 

2.936 X 10'® 

5.0 

0.1140 

6.599 X 10‘® 

0.1141 

6.595 X 10‘® 

10 

0.1621 

9.394 X 10’® 

0.1622 

9.383 X 10'® 

25 

0.2602 

1.515 X 10'^ 

0.2607 

1.511 X 10"^ 

50 

0. 3777 

2.216 X 10“^ 

0.3796 

2.205 X 10"^ 

100 

0. 5652 

3.368 X 10’^ 

0.5720 

3.339 X 10"^ 


The results tabulated in Table 1-2 show that the rms errors of the 
LST with sample-and-hold and the analog controller are very close to those 
of the LST with the controller implemented digitally. 

It is interesting to show that the rms attitude errors of the digital 
LST due to quantization are related to those of the continuous-data system 
in Fig. 1-12. 

Applying the limit as T approaches zero to Eq. (1-31), we get 


ton air, = tim 
T-m T-O 


Mp(z)Mp(z'bz'^dz$^ 


1/2 


(1-46) 


tim M,»(2) = Mf.(s) 
T-vO 


Since 


(1-47) 
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where M|j(s) is given in Eq. (1-23), and Mjj(z) is given by either Eq. 

(1-40) or (1-43), depending on the way the controller is implemented, 
and 

z 'dz = T ds (1-48) 

Eq. (1-46) is written 



r 

[ T 


tcm a*- = 
T^ 

I 1 

2ttj _ 

MD(s)MD(-s)ds«l>D 

-jco 


1/2 


^BD 


(1-49) 


or 

( 1 - 50 ) 

The meaning of this relation is that the mean«square value of the 

attitude error of the continuous -data system is equal to 1/T times the 

attitude error of the digital system as T approaches zero. Table 1-3 

gives the values of and cr|j/ /r<J>.p for various values of T. 

These values are very close to the values of a„r, = 1-605 and = 

^ uU BT 

0.926 X 10 , respectively, especially at very small sampling periods. 

These relations show that, in general, for small sampling periods, 
the attitude error due to a white noise input in a digital system will 
be less than that of the same system without sampling. For example, 
for T = 25 msec, the rms attitude error is 0.2602i/i^ , whereas 
OgD is 1.605v^ . Therefore, the use of the continuous-data LST model 
of Fig. 1-12, as in [1], for the error analysis of the digital LST 
results in conservative results. 


Table 1-3. 


T 

(msec) 

^BD 

/tt„ 

^BT 

0.01 

1.644 

0.927 X 10"^ 

0.1 

1.610 

0.927 X 10"^ 

1.0 

1 .606 

0.928 X 10"^ 

5.0 

1 .612 

0.933 X 10'® 

10 

1.621 

0.939 X 10'® 

25 

1.645 

0.958 X 10'® 

50 

1.689 

0.991 X 10‘® 

100 

1.787 

1.065 X 10'® 


The results in Tables 1-2 and 1-3 again show that for the same 
quantization levels, the quantizer in the displacement channel produces 
far greater attitude errors than those due to the quantizers in the rate 
and torque channels. 



1-6. Effects of Quantization on Po i nting Stability of the Digital Low- 
Cost 1ST - Ouasil inear Analysis, Statistical 
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In this section the effects of quantization on the low-cost LSI 

are studied by means of a quasil inearized equivalent-gain approach [4], 

[5]. Although the analysis is conducted in the statistical sense, there 

is a basic difference between the present analysis and the one conducted 

in Section 1-5. In the previous section, the rms attitude error of the 

LSI is evaluated by treating the quantizers as noise sources that are 

stationary Gaussian processes with zero means. In this section, 

the quantizer will be treated as a nonlinear element whose input is 

a random stochastic process. An equivalent gain, derived for 

the quantizer. The attitude of the LST is then determined with the quantizer 

replaced by K (z), and when the system is subject to a stochastic input, 
eq 

The only restriction with this method is that only one quantizer can be 
considered at a time. Since it has been established that the attitude 
of the LST is more sensitive to the quantizer Qp, we shall consider only 
the system model shown in Fig. 1-14. 

The quantizer Qp is isolated as shown by the block diagram in 
Fig. l-15a; the input is x{t) and the output is y(t). Figure l-15b 
shows the equivalent gain representation of Qp. The equivalent gain 
of Qp is defined as 



(z) 




( 1 - 51 ) 


where 4* (z) denotes the z-transform of the autocorrelation function of 

X X 

x(t), and 4 (z) is the z-transform of the crosscorrelation function of 
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x(t) and y(t). In order to derive K^,q(z). we consider that the input 
to Qg is a stationary Gaussian process with zero mean so that the 
probability function is 


p(x) = 


exp(-x^/2o^) 

ym o 

X 


(1-52) 


where is ttre standard deviation of x(t). 

The crosscorrelation function of x(t) and y(t) is given by 


= j Qp(x)p(x)pxdx 

where 

OpCxCt)] = y(t) 

is the mathematical description of the quantizer, and 

P - 2 


(1-53) 


(1-54) 


(1-55) 


Thus, 




Qp(X)p(X)XdX 


(1-56) 


Taking the z-transform on both sides of Eq. (1-56) and rearranging, we have 

(1-57) 


1 

"eq^^) = ^ = T 

^ XX ' ' 


Qp(X)p(X)XdX 
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Although the equivalent gain is indicated as a function of z, 

the nature of the right-hand side of Eq. (1-57) implies that it is 
always a constant. It is apparent that although our interest is centered 
on the quantizer nonlinearity, in general, the definition of Kgq(^) 
be applied to any coinman nonlinearities found in control systems. 

With reference to the quantizer characteristics of Fig. 1-7, 


o 

II 

-h/2 < X < h/2 

y(t) = -h 

-3h/2 < X £ -h/2 

y(t) = h 

h/2 5 X < 3h/2 

y(t) = -(M-l)h 


y(t) = (N-l)h 

2t!hllh , , < 


. (2N-l)h 

y(t) = -Nh 

-00 < X <^ - ' 2 * " 

y(t) = Nh 

lltMJh < X < » 


where N is a positive integer. 
Equation (1-57) gives 

, fr(2N-l)h/2 

X 


(-Nh)p(X)XdX 


f-(2N-3)h/2 
-(2N-l)h/2 
3h/2 


-(N-Dhp(X)XdX + .. . ^ 
r(2N-l)h/2 


-h/2 


-hp(X)XdX 


h/2 


hp(X)XdX 


(2N-3)h/2 


-3h/2 
(N-Dhp(X)XdX 


(2N-l)h/2 


Nhp(X)XdX 


(1-58) 
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Or, 


1 


X ^ 


+ 2Nh 


2h 


3h/2 

h/2 


p{X)AdA + • • • *■ 2(N-1 )h 


(2N-l)h/2 


p(A)AdA 


1 


(2N-l)h/2 

(2N-3)h/2 


p(A)AdA 


(1-59) 


Since all the integrals ■'n the last equation have the same integrand, 
we evaluate one of the integrals as follows: 


r3h/2 

h/2 


p(\)XdA = 1 

i h/2 


,3h/2 exp{-A^/2o^) 


/2n o. 


AdX 


(1-60) 


Let 


u = X^/2o^ 


Then 


du = 


XdX 


or AdX = du 


and 


3h/2 

h/2 


p(A)XdA 


/2it 


2 2 
9hV8o^ 

2 exp(-u)aj^du 

h /8a^ 


1 




/2it ^ 


exp(- + exp{- -^) 
8»x 


2 1 


(1-61) 


Thus, the equivalent gain for the quantizer with 2N levels of quantization 
is 
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exp(- * exp(- - i) + exp{- * ■■■ * exp(- ) 


-y; ^ 0/ ’ '-''KV 7 ^ ^ ^ 2 

80^ 8,; S 4 80^ 
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Although th8 equivdlent gam has boon dGrived from Qp, the 

response of the closed-loop system of Fig, 1-14 cannot be determined 

by using l^gq(^) directly, since o^, the standard deviation, of x(t) 

is not known. The analysis procedure is outlined as follows. 

For a aiven o « which is the rms value of the input, 

^ c 

(1) The equivalent gain K (z) is computed using Eq. (1-62) for 

various values of a^. 

(2) The values of that correspond to the various values of 
obtained in step (1) are calculated from 






0-63) 


where 


M,(.) 


1 ^ 


0-64) 


and 

\c<"> ' °c 

It is assumed that the input of the system is a white noise. 

The solution of for the given is determined when there is 

a match between from Eq. (1-63) and that used in step (1). 

(3) The rms attitude response of the system is determined from 
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On = 


Ztij 


(j) M(z)M(z‘^ )$j.^(z)z‘^dz 


1/2 


(1-65) 


using the K (z) which corresponds to the obtained in 
eg ^ 

step (2). 

For the 1ST system, a quantization level of h = 0.003 is selected 
for Qp. The quantizer is assumed to saturate after 5000 increments; that 
is, N = 5000 in Eq. (1-62). A digital computer program is prepared 
which automatically cuts off the series of Eq. (1-62) when an additional 
term is contributing less than 10 ^ to the entire result. The sampling 
period T is chosen to be 25 msec, and the system transfer functions 6^(z)» 
Gq(z), and 6 (z) are given by Eqs. (1-38), (1-36), and (1-37), respectively. 
Following the procedure outlined above, the results of the analysis with 
several values of are tabulated in Table 1-4. 

The results in Table 1-4 show that for a given set of quantization 
and saturation levels, both small and large input signals cause the 
quantizer to act as an attenuator. The true characteristics of the system 
as a function of the input are displayed by normalizing Og. The last 
column in Table 1-4 represents the normalized output Ogj^, which is 
defined as 


_ 

^BN - og^o^ 

where 


(1-66) 


‘^Bl " ^^c = ^ 


(1-67) 
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Input 

0 

X 

K 

eq 

No. of terms 
used in K^q(z) 

Output 

No rma 1 1 zed 
output 

^BN 

5x 10"'^ 

5.318X 10"^ 

0.075 

2 

5.20x 10"^ 

0.4 

10'^ 

1.086 X 10"^ 

0.797 

3 

2.40 X 10"^ 

0.92 

10-2 

1 .093 X 10'2 

1.0 

19 

2.68 X 10"^ 

1.0 

10"^ 

1.093X 10’^ 

1.0 

176 

2.60x 10"2 

1.0 

1 

1.093 

1.0 

1753 

2.60 X IQ-^ 

1.0 

10 

10.87 

0.836 

5000 

2.383 

0.916 

50 

53.37 

0.22 

5000 

6.88 

0.529 

Figure 1-16 shows the plot of 

K (z) versus a . 

The plot has the 


significance of a "statistical describing function" of the Quantizer. This 
plot shows that for small inputs, when the input magnitude is comparable to 
the Quantization level, h, the gain drops below unity, and the Quantizer is 
attenuating the signal. For larger inputs the quantizer appears as a unity- 
gain element in a statistical sense, and for very large inputs, where the 
input magnitude is comparable to maximum output (N levels of h), the 
quantizer gain again reduces, and the input is attenuated. 

Figure 1-17 gives a plot of versus which represents the solution 
of Eq. (1-63) and step (1) for the range of considered. 

Figure 1-18 shows a plot of versus o^. It is seen that the normalized 
output reduces from unity in the ranges of when Kg^(z) is less than 1. 

As an interesting comparison, the statistical method of Section 1-5 
assumes Kg^(z) = 1 for all values of and, consequently, it always yields 



Figure 1-16 
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a result of o„., =1. In this sense that method is more conservative and does 

BN 

not TGCognizo the severity of the nonlinedrity due to the quantizer at 
very small and very large input magnitudes. 
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